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The empirical Bath's law states that the average differ- 
ence in magnitude between a mainshock and its largest af- 
tershock is 1.2, regardless of the mainshock magnitude. Fol- 
lowing Vere- Jones' [1969] and Console et al. [2003] , we show 
that the origin of Bath's law is to be found in the selection 
procedure used to define mainshocks and aftershocks rather 
than in any difference in the mechanisms controlling the 
magnitude of the mainshock and of the aftershocks. We use 
the ETAS model of seismicity, which provides a more realis- 
tic model of aftershocks, based on (i) a universal Gutenberg- 
Richter (GR) law for all earthquakes, and on (ii) the increase 
of the number of aftershocks with the mainshock magnitude. 
Using numerical simulations of the ETAS model, we show 
that this model is in good agreement with Bath's law in a 
certain range of the model parameters. 



1. Introduction 

Bath' s law [Bath, 1965] predicts that the average magni- 
tude difference Am between a mainshock and its largest af- 
tershock is 1.2, independently of the mainshock magnitude. 
Many studies have validated Bath's law, with however large 
fluctuations of Am between and 3 from one sequence to 
another one [e.g. Felzer et al, 2002; Console et al., 2003]. In 
addition to providing useful information for understanding 
earthquake processes, Bath's law is also important from a 
societal view point as it gives a prediction of the expected 
size of the potentially most destructive aftershock that fol- 
lows a mainshock. 



2. Vere- Jones' interpretation of Bath's law 

Bath's law is often interpreted as an evidence that main- 
shocks are physically different from other earthquakes and 
have a different magnitude distribution [e.g. Utsu, 1969]. 
In contrast, Vere-Jones [1969] offered a statistical interpre- 
tation, elegant in its simplicity, which consisted in viewing 
the magnitudes of the mainshock and largest aftershock as 
the first and second largest values of a set of independent 
identically distributed (lid) random variables distributed ac- 
cording to the same GR distribution P(m) ~ 10~*™. If the 
same minimum threshold mo applies for both aftershocks 
and mainshocks, this model predicts that Am has the same 
density distribution PAm(Am) ~ ^Q-bAm ^ ^j^^ qj^ distri- 
bution of the sample [Vere-Jones, 1969] with a mean (Am) 
equal to l/(felnlO) ~ 0.43 for b ~ 1. Thus, rather than a 
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distribution peaked at Am « 1.2, Vere-Jones' interpretation 
predicts an exponential distribution with an average signif- 
icantly smaller than Bath's law value ~ 1.2. Such discrep- 
ancies have been ascribed to different magnitude thresholds 
chosen for the definition of mainshocks and largest after- 
shocks and to finite catalog size effects [Vere-Jones, 1969; 
Console et al., 2003]. Improved implementation of Vere- 
Jones' model by Console et al. [2003], taking into account 
the fact that the minimum threshold for aftershock magni- 
tudes is smaller than for mainshocks, has shown that this 
model provides a much better fit to the data, but that 
there is still a minor discrepancy between this model and 
the observations. The results of [Console et al., 2003] for 
a worldwide catalog and for a catalog of seismicity of New 
Zealand are not completely explained by this model, the 
observed value of (Am) being still a little larger than pre- 
dicted. Console et al. [2003] interpret this result as possibly 
due to "a change in the physical environment before and af- 
ter large earthquakes" but they do not rule out the existence 
of a possible bias that may explain the discrepancy between 
their model and the observations. We propose in section 3 
a simple statistical interpretation of Bath's law, which can 
explain this discrepancy without invoking any difference in 
the mechanisms controlling the magnitude of the mainshock 
and of the aftershocks. 

Notwithstanding the appealing simplicity of Vere-Jones' 
interpretation and its success to fit the data, this model does 
not provide a realistic model of aftershocks, and misses some 
important properties of seismicity. In particular, it does not 
take into account the fact that aftershocks represent only a 
subset of the whole seismicity, which are selected as events 
that occurred within a space-time window around and af- 
ter a larger event, called the mainshock, which is supposed 
to have triggered these earthquakes. We first consider as 
a mainshock only the largest earthquake of a catalog of A'^ 
events which have independent magnitudes drawn according 
to the GR law I0~''(™~'"o) ^jth a minimum magnitude mo- 
Only a small subset of size A''aft of the whole catalog occur 
in the specified space-time window used for aftershock selec- 
tion. The largest event in the whole catalog has an average 
magnitude given by {niM) ~ mo -I- (1/6) logjo Let us sort 
the magnitudes of all events in the catalog by descending 
order: mi > m2 > ... > mjv. The largest aftershock, within 
the subset of aftershocks of size A''aft, has an expected over- 
all rank equal to ~ A/A^aft . Using the distribution of the 



magnitude difference mi 



between the largest earth- 



quake (with rank equal to 1) and the event of rank j given 
by [Vere-Jones, 1969] and assuming N ^ Aatt 2> 1, the av- 
erage magnitude difference between the mainshock and its 
largest aftershock is thus given by 



(Am) = (mM - ruA) ~ j logi(,(A/A"aft) (1) 

This expression (1) shows that if the mainshock is taken 
to be the largest event in the catalogue, then the mag- 
nitude difference (Am) is likely to be substantially larger 
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than that predicted by VcTe- J ones' initial formulation. That 
is because the subset formed by the aftershocks is signifi- 
cantly smaller than the original catalogue of size A'^, thus 
the rank of the largest aftershock is not 2 in general. lu 
this sense, the mainshock in a sequence is not a member 
of the set of aftershocks [Utsu, 1969; Evison and Rhoades, 
2001]. The mainshock appears to be an outlier when we 
compare the mainshock magnitude with the aftershock mag- 
nitude distribution, even if all events in the initial cata- 
log have the same magnitude distribution. The fact that 
the mainshock does not belong to the subset of aftershocks 
does not however imply that mainsliocks are physically dif- 
ferent from other earthquakes, in contradiction with pre- 
vious claims of Utsu [1969], but simply results from the 
rules of aftershock selection. Expression (1) retrieves Bath's 
law only for a specific value of the number of aftershocks 

N^H = lQ-H^m)-^Qb(mM-mo) ^^j^Jj ^^^^ ^ 12. 

We use in section 3 the ETAS model of seismicity, in 
order to take into account the rules of aftershock selection 
in time, space and magnitude, and to take into account the 
increase of aftershock productivity with the maisnhock mag- 
nitude. We will generalize the result (1) in the case where 
the mainshock is not the largest event of the whole catalog. 

Using an approach similar to expression (1) and taking 
Bath's law as given led Michael and Jones [1998] and Felzer 
et al. [2002] to deduce that the number of earthquakes trig- 
gered by an earthquake of magnitude m is proportional to 
~ 10™"", with a = b. We shall see below using numeri- 
cal simulations that the ETAS model is also consistent with 
a<b. 

3. Bath's law and the ETAS model 

In order to shed light on the explanation of Bath's law, 
and to investigate the effects of the selection procedure for 
aftershocks, we need a complete model of seismicity, which 
describes the distribution of earthquakes in time, space and 
magnitude, and which incorporates realistic aftershock prop- 
erties. We thus study the Epidemic Type Aftershock Se- 
quence model (ETAS) of seismicity, introduced by [Kagan 
and Knopoff, 1981; Ogata, 1988]. The ETAS model assumes 
that each earthquake triggers aftershocks with a rate (pro- 
ductivity law) increasing as p(m) = KIO"'-"^'"^"^ with its 
magnitude m. A crucial assumption of the ETAS model is 
that all earthquakes have the same magnitude distribution, 
given by the GR law, independently of the past seismicity. 
The seismicity results from the sum of an external constant 
average loading rate and from earthquakes triggered by these 
sources in direct lineage or through a cascade of generations. 

It can be shown that the average total number of after- 
shocks (Aaft(niM)) (including the cascade of indirect after- 
shocks) has the same dependence with the mainshock mag- 
nitude 



(Na.ft{mM)) = 



K 

1-n' 



-10' 



'(mM-mo) 



(2) 



as the immber of direct aftershocks pijUM) given above 
[Helmstetter and Sornette, 2002]. n is the branching ratio, 
defined as the average number of directly triggered earth- 
quakes per earthquake, averaged over all magnitudes. Using 
this model, Felzer et al. [2002] have argued that a must be 
equal to 6 in order to obtain an average difference in magni- 
tude (Am) that is independent of the mainshock magnitude. 
This result is in apparent disagreement with the empirical 
observation a w 0.8 < 6 « 1 reported by Helmstetter [2003] 
using a catalog of seismicity for Southern California. The 
analysis of Felzer et al. [2002] neglects the fluctuations of 



the number of aftershocks from one sequence to another one. 
We have however shown recently [Saichev et al., 2003] that 
there are huge fluctuations of the number of aftershocks per 
sequence for the same mainshock magnitude. We show be- 
low, using numerical simulations of the ETAS model, that 
taking into account these fluctuations has important effects 
on the estimation of (Am) and on its dependence with ttim ■ 
The average magnitude (m^) of the largest event in a 
catalog of iVaft aftershocks with magnitudes larger than mo 
distributed according to the GR law is given by [Feller, 1966] 



(tua) = mo 



_ r jVafl(l 

Jo 



dx 



(3) 



61n(10) 

w mo-|-logio(iVaft)/6 foriVaft>l. (4) 

We derive below an approximate expression for (Am) in the 
ETAS model, which neglects the fluctuations of A'aft, i.e. 
which replaces A'aft by its average value (2) in (4). Using 
this approximation, we obtain 



,., b — a , , 1, / K \ 

(Am) « —^{mM -mo) - - logio ) ' 



(5) 



This approximate relation thus predicts an increase of (Am) 
with the mainshock magnitude if a < 6. Expression (5) thus 
predicts that Bath's law is only recovered for a = 6 [Felzer et 
al., 2002]. Using numerical simulations of the ETAS model, 
we find however that the large fluctuations in aftershock 
numbers due to the cascades of triggered events modify sig- 
nificantly the prediction (5). Adding the constraint that 
aftershocks are usually chosen to be smaller than the main- 
shock further alters the prediction (5). 

We have generated synthetic catalogs with the ETAS 
model to measure (Am) for different values of the main- 
shock magnitude. In this first test, we starts the simula- 
tion with a mainshock of magnitude mm, which generates a 
cascade of direct and indirect aftershocks. We select as "af- 
tershocks" all earthquakes triggered directly or indirectly by 
the mainshock, without any constraint in the time, location, 
or magnitude of these events. For a = 0.8 and 6=1, we 
find that (Am) is much larger than predicted by (5), and 
increases slower with the mainshock magnitude, in better 
agreement with Bath's law than the analytical solution (5) 
(Figure 1). We have checked that the average rmmber of 
aftershocks (A^aft (mM)) is in good agreement with the ana- 
lytical solution (2), and thus a discrepancy with (2) is not 
an explanation for the difference between the results of the 
numerical simulations and the analytical prediction (5). 

The large fluctuations of the total number of aftershocks 
are at the origin of the discrepancy between the observed 
(Am) and the prediction (5), which neglects the fluctua- 
tions of the number of aftershocks. Saichev et al. [2003] 
have recently demonstrated that the total number of after- 
shocks in the ETAS model in the regime a > b/2 has an 
asymptotic power-law distribution in the tail with an expo- 
nent of the cumulative distribution smaller than 1, even in 
the subcritical regime (defined by a branching ratio n < 1). 
These huge fluctuations arise from the cascades of trigger- 
ing and from the power-law distribution of the number of 
triggered earthquakes per triggering earthquake appearing 
as a combination of the GR law and of the productivity 
law p(m). Practically, this means that the aftershock num- 
ber fluctuates widely from realization to realization and the 
average will be controlled by a few sequences that happen 
to have an unusually large number of aftershocks. Numer- 
ical simulations show that (Am) is not controlled by the 
average number of aftershocks, but by its "typical" value, 
which is much smaller than the average value. Therefore, 
the expression (5) of (Am) obtained by replacing TVaft by its 



HELMSTETTER AND SORNETTE: BaTH'S LAW, GUTENBERG-RICHTER LAW AND AFTERSHOCK PROPEmES 



1.6 
1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 


-0.2, 
-0.4 



AM=1.2 _j 

n O 
O ° ° ^ 

^ooooooooO° 

O ^-"^ 
o ^--^ 




- 


o 




o 













0123456789 
mainshock magnitude 

Figure 1. Average magnitude difference (Am) between 
the mainshock and its largest aftershock as a function 
of the mainshock magnitude (open circle), for numerical 
simulations of the ETAS model with parameters n = 0.8, 
a = 0.8, 6 = 1, mo — and with an Omori exponent 
p = 1.2. The continuous line is the prediction using the 
approximate analytical solution (5) for (Am). 



average value (2) in (3) is a very bad approximation. Using 
the exact distribution of the number of aftershocks given in 
[Saichev et al, 2003], we can obtain the asymptotic expres- 
sion for large itim, which recovers the dependence of (Am) 
with ruM predicted by (5). 

For large mainshock magnitudes, the relative fluctuations 
of the total number of aftershocks per mainshock are weaker. 
Therefore, the obtained average magnitude difference (Am) 
tends to recover the linear dependence (5) with the main- 
shock magnitude, represented by the continuous line in Fig- 
ure 1. Our numerical simulations show that a constant value 
of (Am) in a wide range of magnitudes can be reproduced 
using the ETAS model if a < 6. Our results also predict 
that Bath's law should fail for large mainshock magnitudes 
according to (5) if a is smaller than b. It is however doubtful 
that this deviation from Bath's law can be observed in real 
data as the number of large mainshocks is small. 

While the impact of fluctuations in the number of after- 
shocks produces a value of (Am) larger than predicted by 
(5) and in better agreement with Bath's law, the average 
magnitude difference (Am) « 0.7 remains smaller than the 
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Figure 2. Average magnitude difference (Am) between 
a mainshock and its largest aftershock, for numerical 
simulations of the ETAS model with parameters 6=1, 
c = 0.001 day, p — 1.2, a minimum magnitude mo = 2, 
a maximum magnitude rrimax = 8.5 and a constant load- 
ing /i = 300 events per day. Each curve corresponds to 
a different value of the ETAS parameters: a = 0.8 and 
n = 0.76 (crosses), a — 0.5 and n = 0.8 (diamonds) and 
a = 1 and n = 0.6 (circles). The error bars gives the un- 
certainty of (Am) (1 standard deviation). The horizontal 
dashed line is the empirical value (Am) = 1.2. 



empirical value (Am) ~ 1.2. However, we have not yet taken 
into account the constraints of aftershocks selection, which 
will further modify (Am). In the simulations giving Figure 
1, all earthquakes triggered (directly or indirectly) by the 
mainshock have been considered as aftershocks even if they 
were larger than the mainshock. In real data, the difficulty 
of identifying aftershocks and the usual constraint that af- 
tershocks are smaller than the mainshock can be expected to 
affect the relation between (Am) and the mainshock mag- 
nitude. The selection of aftershocks requires the choice of 
a space-time window to distinguish aftershocks from back- 
ground events. A significant fraction of aftershocks can thus 
be missed. As a consequence, the value of Am will increase. 

In order to quantify the impact of these constraints, we 
have generated synthetic catalogs using the ETAS model, 
which include a realistic spatio-temporal distribution of af- 
tershocks. Specifically, according to the ETAS model, the 
number of aftershocks triggered directly by an event of mag- 
nitude m, at a time t after the mainshock and at a distance 
r is given by 



0™(f,r) = ni^-^lO"f'"-'""' 
6 



where n is the branching ratio, p is the exponent of the lo- 
cal Omori's law (which is generally larger than the observed 
Omori exponent) and dm is the characteristic size of the 
aftershock cluster of a magnitude m earthquake given by 
dm = 0.01 X lO"-^" km. 

We have then applied standard rules for the selection of 
aftershocks. We consider as a potential mainshock each 
earthquake that has not been preceded by a larger earth- 
quake in a space-time window Rc x Tc- This rule allows us to 
remove the influence of previous earthquakes and to obtain 
an estimate of the rate of seismicity triggered by this main- 
shock. The constant i?c is flxed equal to the size ~ 100 km 
of the largest cluster in the catalog and Tc = 100 days. We 
then define aftershocks as all events occurring in a space 
time window R{mM) x T{mM) after a mainshock of magni- 
tude niM, where both R{mM) = 2.5 x IO'^'^^a^ km and 
T{mM) = 10/3 X lo(2/3){'"Af -5) (jayg increase with the main- 
shock magnitude ttim [Kagan, 1996; Console et al, 2003]. 

The results for different values of a are represented in Fig- 
ure 2. For intermediate mainshock magnitude, the average 
magnitude difference (Am) for a = 0.8 is significantly larger 
than found in Figure 1 without the selection procedure, be- 
cause mainshocks which trigger a larger event are rejected, 
and because the rules of selection (with a time-space window 
R{m) and T(m) increasing with m) reject a large number of 
aftershocks, especially for small mainshocks. For small mag- 
nitude niM, (Am) is small and then increases rapidly with 
m. This regime is not pertinent because most mainshocks 
do not trigger any aftershock and are thus rejected from 
the analysis. Most studies have considered only mainshocks 
with magnitude m > mo + 2, where mo is the minimum 
detection threshold. For a = 0.8 or a = 1, the magnitude 
difference is ~ 1.2 in a large range of mainshock magnitudes, 
in agreement with Bath's law. For a = 1, there is a slight 
decrease of (Am) with ium- For a = 0.5, we observe a fast 
increase of (Am) with niM, which is not consistent with the 
observations of Bath's law [e.g. Felzer et al., 2002; Console 
et al., 2003]. The shape of the curves (Am) is mostly con- 
trolled by a. The other parameters of the ETAS model and 
the rules of aftershock selection increase or decrease (Am) 
but do not change the scaling of (Am) with the mainshock 
magnitude. 

4. Discussion and conclusion 

We have first shown that the standard interpretation of 
Bath's law in terms of the two largest events of a self-similar 
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set of independent events is incorrect. We have stressed the 
importance of the selection process of aftershocks, which 
represent only a subset of the whole scisrnicity catalog. Our 
point is that the average magnitude difference (Am) is not 
only controlled by the magnitude distribution but also by 
the aftershock productivity. A large magnitude difference 
(Am) can be explained by a low aftershock productivity. 

Using numerical simulations of the ETAS model, we have 
shown that this model is in good agreement with Bath's law 
in a certain range of the model parameters. We have pointed 
out the importance of the selection process of aftershocks, of 
the constraint that aftershocks are smaller than the main- 
shock and of the fluctuation of the number of aftershocks per 
sequence in the determination of the value of (Am), and in 
its apparent independence as a function of the mainshock 
magnitude. In the ETAS model, the cascades of multiple 
triggering induce large fluctuations of the total number of 
aftershocks. These fluctuations in turn induce a modifica- 
tion of the scaling of (Am) with the mainshock magnitude 
by comparison with the predictions neglecting these fluctu- 
ations. The constraints due to aftershock selection further 
affect the value of (Am). Observations that (Am) does not 
vary significantly with the mainshock magnitude requires 
that the exponent of the aftershock productivity law is in 
the range 0.8 < a < 1. Bath's law is thus consistent with the 
regime a < 6 in which earthquake triggering is dominated 
by the smallest earthquakes [Helmstetter, 2003]. 
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